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Abstract — Stopping sets, and in particular their numbers and 
sizes, play an important role in determining the performance of 
iterative decoders of linear codes over binary erasure channels. In 
the 2004 Shannon Lecture, McEliece presented an expression for 
the number of stopping sets of size three for a full-rank parity- 
check matrix of the Hamming code. In this correspondence, we 
derive an expression for the number of stopping sets of any given 
size for the same parity-check matrix. 

Index Terms- Hamming code, linear code, parity-check matrix, 
stopping set enumerator, weight enumerator. 

I. Introduction 

Let C be a linear binary [n, k, d] block code, where n, k, and 
d denote the code's length, dimension, and Hamming distance, 
respectively. The set of codewords of C can be defined as the 
null space of the row space of an r x n binary parity-check 
matrix H — {hij) of rank n — k. Assuming all rows in H are 
different, n - k < r < 2""'=. 

Let iS be a subset of {l,2,...,n} and T be a subset 
of {l,2,...,r}. For any H = (hij) of size r x n, let 
(hij) where i e T and j G S. Then, Hj is a |T| x |iS| 



submatrix of H. For simplicity, we write and H to 
denote Hj in case T = {1, 2, . . . , ?■} and S — {1,2,..., n}, 
respectively. 

A set S is the support of a codeword if and only if all rows 
in Hg have even weight, i.e., if and only if 



\{jeS:h,,, = l}\=0{2) 



Vi = 1,2,. 



(1) 



A set iS is a stopping set if and only if does not contain 
a row of weight one, i.e., if and only if 



|{je5:/i,,, = 1}|^1 



V* = l,2,., 



(2) 



The polynomial A{x) — J2?=o ^i^;', where Ai is the num- 
ber of codewords of weight /, is called the weight enumerator 
of code C. It holds that 

d = min{/ > 1 : A; > 0}. 
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The polynomial S{x) = X^iLo ^ix^, where Si is the number 
of stopping sets of size I, is called the stopping set enumerator 
of parity-check matrix H. Let s denote the smallest size of a 
non-empty stopping set, i.e., 

s = min{/ > 1 : 5; > 0}. 

Notice from ([l) and (|2j that the support of any codeword 
is a stopping set. Therefore, Si > Ai for I = 0,1, . . . ,n and 
s<d. 

Considering the vacuous case in which I = 0, we notice, 
from the definitions, that the empty set is both the support of 
a codeword and a stopping set for any code and any parity- 
check matrix. Hence, Aq = Sq = 1. Furthermore, from the 
observation that ([0 and (|2j are equivalent for sets S with 
|5| < 2, it follows that Ai = Si for / < 2. In particular, 
5o = 1 and 5i = 52 = for any parity-check matrix of a 
code of minimum distance d > 3. On the other hand, if every 
row in H has weight at least equal to ri — / + 2 and 5 is a 
subset of {1,2,..., n} of size I, then every row in Hg has 
weight at least equal to two. Hence, 5 is a stopping set of H 
and 

Si^i':). (3) 



n\ 
l) ' 

The notion of stopping sets is important in the context of 
iterative decoders (using H-based Tanner graphs), in particular 
for low-density parity-check codes [3]. For example, on the 
binary erasure channel, an iterative decoder will not lead to 
successful decoding if and only if the set of erased positions 
contains a non-empty stopping set. Hence, the minimum (non- 
empty) stopping set size s and the cardinality Sg are perfor- 
mance indicators for iterative decoding, alike the minimum 
distance d and the number of minimum weight codewords 
for maximum-likelihood decoding. However, contrary to the 
weight enumerator, which is fixed for a code C, the stopping 
set enumerator depends on the choice of the parity-check 
matrix H. 

In the 2004 Shannon lecture, McEliece [6] presented the 
following expression for the number of stopping sets of size 



three in a [2™ - 1,2^ 



1, 3] Hamming code: 



^3 



-)m+l\ 



1 



2.322 



(4) 



He did not mention expHcitly which parity-check matrix he 
had in mind, but from the context, it was clear that it was the 
full rank m x (2™ — 1) parity-check matrix. This is the parity- 
check matrix of minimum number of rows for the Hamming 
code. The value of S3 is of particular interest since 5*0 = 1 and 
S*! = 5*2 = as the Hamming code has minimum distance 
three which implies that s = 3. Hence, the value of ^3 given 
in corresponding to the full rank parity-check matrix of the 
Hamming code can be used to obtain a good estimate of the 
performance of the simplest iterative decoder for the Hamming 
code. However, as the multiplicities Ai of codewords in a 
code of weight I > d influence the performance of maximum- 
likelihood decoding, the multiplicities Si of stopping sets of 
size I > s also influence the performance of iterative decoding. 
In this correspondence, we derive the values Si for all I for 
the full-rank parity-check matrix of the Hamming code. 

Section II gives an expression for Si in terms of certain 
parameters that can be derived from the parity-check matrix of 
any linear code. In Section III, we determine these parameters 
for the full-rank parity check matrices of Hamming codes and 
derive expressions for their stopping set enumerators. 

II. Computing the Stopping Set Enumerator for 
Linear Codes 

In this section, we derive an expression for the coefficients 
of S{x) for an r X n parity-check matrix H of an [n, k, d] 
binary linear code C. 

For each subset T of {1,2,..., r}, we define zj- to be the 
number of columns in of weight zero. In particular, if 
T is empty, then zq- = n. If 3^ is a subset of {1, 2, . . . , n} 
such that 1 3^ I = p, each row in Hy has weight one, and each 
column in Hy has non-zero weight, then we say that y is 
of type p with respect to T. Let Y{T,p) be the number of 
subsets of {1, 2, . . . , n} of type p with respect to T. Since 
every column in Hy should have non-zero weight and every 
row should have weight one, it follows that 



YiT,p) 



ifr = 0Ap>lorr7^0Ap = O, 

1 ifr = 0Ap = O. 



(5) 

Example 1: Consider the full-rank parity-check matrix H 
of the [7, 4, 3] Hamming code given by 



H 



10 10 10 1 
110 110 
11110 



Let T = {2,3}. Then, zr = 1 as has exactly one all-zero 
column and Y{T, 2) = 4 since there are exactly four subsets 
of type p — 2 with respect to T = {2, 3} given by {3, 5}, 
{3,6}, {4,5}, and {4,6}. Indeed, is a p x \T\ matrix 
each row of which has weight one and each column of which 
is nonzero if and only if y is one of these four subsets. 



Theorem 1: Let H be an r x n parity-check matrix of a 
linear code C of length n. Then, for Z = 0, 1, . . . , n. 



Si = 



rc{l,2,...,r} p=0 ^ ^ 

Proof. It follows from Q that the result holds for I = Q. 
Therefore, we assume I >\. Let 5 be a non-empty subset of 
{1,2,..., n}. Following [8], we say that row h in H covers S 
if has weight one. For i — 1, 2, . . . , r and I — 1,2, ... ^n, 
let Q\ and Q\ be the family of subsets of {1, 2, . . . , n} and 
size I that are covered and not covered, respectively, by the ith 
row in H. A subset iS C {1, 2, . . . , n} of size I is a stopping 
set if and only if it is not covered by any row of H, i.e., 
S e n^^iOf. In particular. Si = | n[^i Q[|. The Principle of 
Inclusion and Exclusion (see, e.g., [2], Theorem B, p. 178), 
yields 



i\ - 



rc{i,2,.. 



(-i)i^'ln 



(6) 



In the above sum, T runs over all 2^ subsets of {1, 2, . . . , r} 
where, for the term corresponding to the empty set 0, 



n.ee Q\\ = 



which is the total number of subsets of {1,2,..., n} of size I. 
Notice that a set S of size I belongs to fligrQi if and only if 
each row in has weight one. This is the case if and only 
if, for some p = 0,1, . . . ,1, S contains a subset y of type p 
with respect to T and all the columns indexed by S\y have 
zero weights. Therefore, 

p=0 ^ ^ 

Combining this with (|6j, the proof is complete. ■ 

III. Computing the Stopping Set Enumerator for 
Hamming Codes 

A full-rank parity-check matrix of a [2™ — 1, 2™ — m — 1, 3] 
Hamming code, where m > 3, is an m x (2™ — 1) matrix 
whose columns are the distinct non-zero vectors of length to. 
From now on, we take H to be such a matrix with r = m 
and n — 2™ — 1, and 5*; to be the number of stopping sets of 
H of size I. In particular, since every row in H has weight 
2™~i, we have from (|3} 

for 2™-i + 1 < / < 2™ - 1. Our goal is to determine Si for 
all I. 

In the following derivations, we make use of Stirling num- 
bers. Following the notation of Comtet [2] and Riordan [7], we 
denote by s{n, k) and S{n, k) the Stirling numbers of the first 
kind and of the second kind, respectively. Notice that n and 
k are not necessarily the length and the dimension of a code. 
For n > 1, (— l)"^'"'s(n, fc) is the number of permutations of 
n elements which have exactly k cycles and S{n, k) is the 



2 



number of ways of partitioning a set of n elements into k 
non-empty subsets, see, e.g., [1], p. 824, where (fTl- dlOt can 
also be found. We also define s(0, 0) = S'(0, 0) = 1. Notice 
that for n > 1, both s(n, k) and S{n, k) are equal to zero if 
A: < or fc > n, see e.g., [7]. 

Stirling numbers of the first kind satisfy the following 
polynomial identity in x 



s(n, k)x^ 



n 



fc=o 



(7) 



and the recursion 

s{n + 1, fc) = s{n, fc — 1) — ns{n, k), (8) 



for n > k > 1. We also have 

s(7i,l) = (-ir-i(n-l)!. 



(9) 



Stirling numbers of the second kind can be computed 
expUcitly using 



i=0 



(10) 



The reader can find a lucid treatment of Stirling numbers of 
both kinds in [4] but with different notation and where Stirling 
numbers of the first kind are defined as (— l)"^'''s(n, k). 

Lemma 1: For n > 0, the following polynomial identity in 
X holds 

^ Y,s{n + l,k)x''-\ 



k=l 



Proof. From the definition of binomial coefficients, we have 

' X — 1 \ X — n 



C) 



n / X \n 
Combining this with we get 

a; — 1\ X — n 



n 



■ ^ s{n, k)x^ 



k=0 



k=0 k=0 

n+l n 

= y^s(n, fc- l)x''~^ -ny^s(n, fc)x''"^ 

fe=l k=0 

n+l 

= {s(n, fc - 1) — ns{n, fc)) x''~^ 

k=l 

since ns{n, 0) = and s{n, n + 1) = for n > 0. The result 
now follows from the recurrence relation (|8j. ■ 
Lemma 2: Let H be an m x (2™ — 1) parity-check matrix 
of a Hamming code of length 2'" — 1. Then, for any set T C 
{1,2,..., to} and any p = 0, 1, . . . , n, 

zr = 2™- 1^1 - 1 



and 



y(r,p) = s{\rip) 2("-i^i)p. 



Proof. Each non-zero vector of length to appears exactly once 
as a column in H. Hence, every zero-weight vector of length 
T| appears exactly 2™^!^' — 1 times in H^. This verifies 
the expression for zq-- From (|5} and ( II 0> . it follows that the 
second part of the lemma holds in case T = 0. Therefore, 
assume in the following that T 7^ 0. Hence, 1 < |T| < to. 
For each j G {1, 2, . . . , 2" - 1}, let 



{leT: h,j = 1}, 



i.e., Hj is the set of indices of ones in the jth column of H . 
Notice that a subset y of size p is of type p with respect to T 
if and only if Hy is a |T| x p matrix, each column of which 
has non-zero weight, and each row of which has weight one. 
This is the case if and only if, for j E y, the sets HJ are 
non-empty, disjoint, and their union is T. There are iS'(|T|,p) 
ways to partition a set of |T| elements into p non-empty 
disjoint subsets. Each non-zero vector of length |T| appears 
exactly 2'"~l^l times as a column in H-^. Therefore, for each 
partition, there are 2™~l^l choices for j, corresponding to 
indices of identical non-zero columns in H^, such that TiJ is 
one of the p subsets in the partition. Since there are p subsets 
in each partition, and there are S{\T\,p) partitions, the total 
number of |T| x p submatrices in H^, where < p < I, 
each of its columns has non-zero weight and each of its rows 
has weight one, is S{\T\,p) 2(™-l^l)P. The number of such 
matrices equals Y{T,p). ■ 
Lemma 3: Let H be an to x (2™ — 1) parity-check matrix of 
a Hamming code of length 2™-l. Then, for ? = 0, 1, . . . , 2™- 
1, 



{m — i)p 



1 



i-p 



t=0 p=0 

Proof. From Lemma |2l it follows that for fixed to and p, 
both z-T and Y{T,p) depend only on the cardinality of T. 
Since there are (™) subsets T of cardinality t, then the result 
follows by using the expressions derived for zt and Y{T,p) 
in Theorem [n ■ 
Lemma |3] when combined with (I10> . gives an explicit 
expression for Si for / = 0, 1, . . . , 2™ — 1. Notice that for 
each value of /, Si depends only on to. However, contrary to 
the expression in 0, the expression in the lemma does not 
show clearly this dependency. In the next lemma, we make the 
dependency clear For this purpose, we define for nonnegative 
integers q and v, 

b{q,v) - ^(-1)P s{p +l,p-q + l). (11) 

p=0 

Lemma 4: Let H be an to x (2™ — 1) parity-check matrix of 
a Hamming code of length 2™ — 1. Then, for 1 = 0,1,..., 2"' — 
1, 

= 7i E E(-i)" ( ^ ) ^(9, v){2'-'^ -{I- v)r. 



3 



Proof. From Lemma ^ we have for p < I and t < m 



j(m — (li — 1) 



1 

^ q=0 

2(m-t){l~q-p) 

1 ' 

^ q=0 

_2im-t){l-q-p) ^ (22) 

since s{n, k) ~ for fc < < n. Also, from (I1Q> . we obtain 



-1)' 



E (-1)^' 



l-v 



l-p 



E(-i)'" 



(Z-t;)*. (13) 



smce 



l-v 

for V < I — p. Substituting ( I12t and ( I13> in the expression of 
S"/ given in Lemma |3j we get 



a = E(-i)*(T)i: 



'"^-'-^B-ir(,!„)('-.)' 

ty{m — t)p ^ 

^sil-p+l,l-q-p+ l)2('"-*)('-«-P) 



t=0 p=0 D=0 



/ - V 



• - w)* ^ s(/ - p + 1 J - g - p + i)2(™-*)('-'?) . 
Noticing that 

we obtain by interchanging orders of summations 

- nEB-ir(:)B-i)'-'(,:, 

5=0 v=0 ^ ^ p=0 ^ ^ 

•s(?-p+l,/-g-p + l) 



E(-1)'"M,! )sa-p+l,/-g-p+l). 

p=0 ^ 

Replacing p by I — p and noticing that 



= 



whenever p > u, it follows that the last sum is identical to 
h{q,v). ■ 

To obtain an expression for v) that does not involve 
Stirling numbers, we derive a recursion for v) along with 
boundary values to apply the recursion. 

Lemma 5: For nonnegative numbers q and v, we have 

if q > V or q = < V, 
ql ii q = V, 

and for q > 1 and v > 2, 

b{q, v) = vb{q - l,v - 1) - {v - l)b{q -l,v- 2). 
Proof. For < p < u < g, we have p — q + 1 < and, 
therefore, s{p + l,p — q + 1) = 0. From the definition of 
b{q, v) given in il 1> . it follows that b{q, v) = in case q > v. 
In case g = 0, it follows from the definition of b{q, v) and the 
fact that s(p+ l,p+ 1) = 1 that b{0, v) ~ whenever v > 1. 
In case q — v, we have 

b{q,q)^J2(~ir('^)s{p+l,p-q+l) 

p=0 

^{~irs{q+l,l)=q\ 

using (|9jl and since s{p + l,p — q + 1) = for p — q + 1 < 0. 
Next, we prove the recursion for q > 1 and i; > 2. Using 
the recurrence formula for the binomial coefficients followed 
by the recurrence formula (|8j for Stirling numbers of the first 
kind, we obtain 

' V — l\ f V — V 



b{q,v)=Y.^-lY 



p=0 



p 



■ s{p + l,p - g + 1) 
p=0 \ P / 



E(-i)' 



p=0 



V - 1 
p-1 



• (s(p, p - g) - ps(p, p - g + 1)) 

= -E(-lf (''Zj)p-*(P'P-9+l) (14) 
p=o ^P ^ 

=eWC)(p+i) 



s(p+l,p-(7 + 2) 



■u-l 



5](-l)*( ")2(™-*)('-^)G-^;)* 



u - 1 



= E(-1)M^ ^ ),(p+l,p_5 + 2) 
p=0 ^ P 



?fEE(-irf!)(2'--('-)) 



5=0 D=0 



^(-i)pr ^ )ps(p+i,p-g+2). 

p=0 ^ P 
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TABLE I 
b{q,v) FORO <q,v<7 



q\v 





1 


2 


3 


4 


5 


6 


7 





1 























1 





1 


-1 

















2 








2 


-5 


3 











3 











6 


-26 


35 


-15 





4 














24 


-154 


340 


-315 


5 

















120 


-1044 


3304 


6 




















720 


-8028 


7 























5040 



Next, we use 



followed by (|8} to conclude that 





















v-1 



v-2 
p-l 



v-2 
p-l 



b{q, v) ^b{q -l,v-l) + iv- 1) Y.{-ir 

■s{p+ l,p- q + 2) 

= b{q-l,v-l) + {v-l)Y,{-ir 

■{s{p,p - q + 1) -ps{p,p-q + 2)) 
= b{q-l,v-l) 

-iv-l) (''~'^)s{p+l,p-q + 2) 

-{v-1) 1 1) p'(p^p -9+2) 

= b{q-l,v-l)-{v- l)b{q -l,v- 2) 
+ {v-l)b{q-l,v-l), 

where we used (II 1> and J14> after replacing v hy v — 1 and q 
by g — 1. Hence, we have 

b{q, v) = vb{q -l,v -1) - {v - l)b{q -l,v- 2). 

■ 

The values of b{q, v) for < q,v < 7 are listed in Table I. 

Now we state an explicit expression for b{q,v). For this 
purpose, we write a <2 b for real numbers a and b if and only 
if 6 — a > 2. For example, 3 <2 5 but 3 ■iit2 4. 

Lemma 6: For nonnegative integers q and v > q, we have 



biq,v)^{-l)^-H,\ 



0=fco<2fel<2---<2'£u-g<2fei.-<j+l=l' + 2 



n 

i=l 



(15) 



Proof. Let b{q, v) denote the expression in the right hand side 
of il5\ . It suffices to verify that b{q, v) satisfies the recursion 
given in Lemma |5] for q > 1 and v > 2 and that b{q, v) and 
b{q, v) agree in the cases q ^ < v and q = v. First, we notice 
that in case = q < v, the summation in the expression of 



b{q, v) is over all tuples — ko, ki, . . . , ky, kv+i = v + 2 for 
which ki — > 2 for i — 1,2, . . . ,v + 1. For such tuples 

v+l 

v + 2^ fc„+i -ko= ^{h - fc,_i) > 2{v + 1) > V + 2. 

4=1 

This contradiction shows that the sum in the expression of 
b{q, v) is empty and, therefore, 6(0, f ) = for v > 1, which 
is identical to b{Q,v) given in Lemma|5] Next, we notice that 
in case v = q, the summation in the expression of b{q, v) 
is over one tuple only, (fco, fci) with fco = and fci —q + 2, 
while the product is empty, and hence equals one. This verifies 
that b{q, q) = q\, which is also identical to b{q, q) given in 
Lemma |5] Finally, we verify the recursion 

b{q, v) = vb{q -l,v-l)-{v- l)b{q -l,v- 2) (16) 

for q > 1 and v > 2. We notice that 

vb{q - 1, t; - 1) = v{-lY-i{v - 1)! 



E 



v-q 



0=ko<2ki<2---<2ki,-q<2k^- 



i=v+l i=l 



(-1) 

E 



0=fco<2fel<2---<2fcii-(j<2feu--g + l=f + 1 



and that 



{v - l)b{q -l,v-2) = {v- l){-iy-'i-\v - 2)! 

v-q-l 

^ n ^ 

0=feo<2'£l<2---<2feii-g=f * = 1 

^-{-iy~iv\ 

v-q-l 

S n^.- OS) 

0=/co<2fel<2---<2fcu-n=t) 



We notice that the sets of tuples kQ,...,ky^q over which 
the summations in (I17> and ( I18> run are disjoint and their 
union gives the set of tuples over which the summation in the 
expression of b{q, v) runs. Hence, M6\ holds. ■ 

The following corollary gives values of b{q, v) for selected 
pairs [q, v). 

Corollary 1: For nonnegative integers q. 



^9,9 + 1) 



9+1 , 



k=2 



9-1 



biq,2q) = {-irY[{2q-2i-l) 

1=0 

b{q,v) — for t; > 2g. 
Proof. The statement regarding b{q, q) follows from Lemma|5l 
The expression for b{q,q + 1) follows readily from M5\ . For 
V > 2q, we notice that the summation in the right hand side of 
(I15> runs over all tuples = fco, fci, ... , ky-q, ky-q+i = v + 2 



5 



for which ki 
tuples, 



fci_i > 2 for i = 1,2, 



■,v - q + 



1. For such 



-9+1 



l-q+1 

E 

1=1 



i)>2iv-q+l). 



Hence, for v > 2q, the sum in the expression of b{q, v)is 
empty and b{q, 2v) = 0. For v ~ 2q, the summation is over 
one tuple only for which ki — 2i for i — 0,1, . . . ,q+l. This 
proves the expression of b{q, 2q). ■ 

Combining Lemmas 0]|5l and|6l we obtain the major result 
of this correspondence. 

Theorem 2: Let H be an m x (2™ — 1) parity-check matrix 
of a Hamming code of length n — 2™ — 1. Then, for / = 
0, 1,...,2™-1, 



Si 



I min{2g,/} 
q=0 



v—q 



("ir 



b{q,v){2^-'^ -{i-v)y 



A, 



where b{q,v) is given in Lemmas [S] and |6l 

Corollary 2: Let H be an to x (2™ — 1) parity-check matrix 
of a Hamming code of length n = 2™ — 1. Then, for a fixed 
value of / > 3, Si behaves asymptotically as a function of the 
code length n as 

Si . (19) 

It is interesting to compare the expression of Si with that 
of Ai, the number of codewords of weight / in a Hamming 
code of length n = 2™ — 1, given by 

This expression can be verified by noticing that A^ = \, Ai = 
0, and that it satisifies the recursion stated on page 129 of [5]. 
For a fixed value of Z > 3, Ai behaves asymptotically as a 
function of the code length n as 

The asymptotic expressions of Si and Ai for / > 3 can be 
directly justified as follows. First notice that l\Si is the number 
of TO X ^ binary matrices with different non-zero columns such 
that no row has weight one. The total number of to x I binary 
matrices without any row of weight one is 

(2'-0"- Amongst 

these matrices, there are at most 12^'-^^^'^ matrices with at least 
one all-zero column and at most (2)2'^'^^^™ matrices with at 
least two identical columns for I > 2. Therefore, by the union 
bound we get 



(2' - /)" - 



l)m 



< usi < {2^ -ly 



Hence, for a fixed / > 3, as n = 2™ — 1 tends to infinity, the 
ratio l\Si/(2' - I)"' tends to one, and (O follows. 

A similar argument leads to i20\ . Notice that II Ai is the 
number of to x Z matrices with different non-zero columns and 
even-weight rows. The total number of to x Z binary matrices 
with even-weight rows is 2('^^)™. Amongst these matrices. 



there are at most Z2^'^^)™ matrices with at least one all-zero 
column and at most (2)2('^'^^™ matrices with at least two 
identical columns for I > 3. Therefore, by the union bound 
we get 

V 



< l\A, < 2('-i) 



Hence, for a fixed Z > 3, as n = 2™ — 1 tends to infinity, the 
ratio /!Ai/2('~^'" tends to one, and (|20j follows. 

Direct application of Theorem |2] gives the following results 
including the expression of ^3 in stated by McEliece [6] 



So- 
S3 

84^ ■ 



l,Si = S2 = 0, 

- (5" - 3 x 3™ + 2 X 2"*) . 
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1 

24 



(12'^ 



6x6*^ 



4 X 5" + 3 X 4'= 



1 

120 



-20 X 3'= 
(27™ - 



14 X 2™) 



10 X 13™ - 5 X 12" 

+15 X 7" + 50 X 6" + 20 X 5" 
-35 X 4™ - 130 X 3" + 94 X 2™) 



Furthermore, if Sm{x) denotes the stopping set enumerator 
for the full-rank parity-check matrix of a Hamming code of 
length 2™ - 1, then 

Ssix) ^ 1 + lOx^ + 23x^ + 2lx^ + 7x^ + x\ 
Si{x) = 1 + 69x^ + 526a;^ + 1979a:^ + 43332:^ + 6211x^ 
+6403x8 + 5005x^ + 300Sx^° + 1365a;" 
+455x^2 + 105a;^^ + 15a;^'' + a;^^ 
Ssix) = 1 + 410x3 + 8215x'* + 83590x^ + 519481x^ 
+2243175x^ + 7378485x8 + 19645915x^ 
+43951765x^° + 84432075x" + 141011325x^2 
+206216675x^3 + 265174125x1"* + 300538995x1^ 
+300540115x16 + 265182525x1^ + 206253075x1* 
+ 141120525x1^ + 84672315x2" + 44352165x^1 
+20160075x^2 + 7888725x23 + 2629575x2" 
+736281x25 + 169911x26 + 31465x2^ 



+4495x^'' + 465x^^ + 31x3" _^ ^31^ 
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